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Abstract 

In  this  study,  numerical  analyses  were  performed  to  predict  the  initial  crack  length  of 
specimens  with  two  different  sizes  of  holes  using  multi-level  simulation  techniques. 
The  specimens  were  subjected  to  a  constant  strain  rate  of  0.33  in/in/min  at  room 
temperature.  The  criterion  for  determining  the  initial  crack  length  was  based  on  the 
instability  of  the  material  near  the  edge  of  the  hole.  The  results  of  the  analyses  show 
that  the  predicted  initial  crack  lengths  compared  well  with  the  experimental 
measurements. 

Introduction 

In  this  study,  a  micro-macromechanical  approach  was  used  to  predict  the  initial  crack 
length  near  the  edge  of  the  hole  in  solid  propellant  specimens.  The  approach  was 
based  on  a  simplified  micromechanical  model,  damage  mechanics  at  the  micro-level, 
and  finite  element  analysis  at  the  macro-level  (1-3).  Both  micromechanical  and 
macromechanical  analyses  were  conducted  in  tandem.  The  developed  technique^ 
together  with  a  mechanistic  criterion^ was  used  to  predict  the  initial  crack  length  in 
high  stress  regions.  The  criterion  was  based  on  the  instability  of  the  damaged  material 
just  ahead  of  the  crack  tip.  The  initial  crack  length  is  equal  to  the  length  of  unstable 
material  zone  when  the  damage  at  the  crack  tip  element  is  saturated.  Based  on  the 
definition  of  the  initial  crack  length  and  the  micro-macromechanical  approach,  the 
initial  crack  lengths  in  the  high  stress  regions  were  predicted.  The  predicted  initial 
crack  lengths  and  the  experimentally  measured  values  were  compared  and  the  results 
are  discussed. 
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Numerical  Modeling 


The  computer  modeling  and  simulation  technique  used  in  this  study  is  called  a 
micro/macro-approach.  The  approach  is  similar  to  the  local-global  approach  in 
structural  analysis.  The  micro/macro-approach  utilizes  the  two  levels  of  analyses: 
micro-analysis  and  macro-analysis.  The  micro-analysis  is  performed  at  the  constituent 
material  levej^hkg)  the  reinforcing  (particle)  material  and  the  binding  matrix  materialA 
while  the  macro-analysis  is  undertaken  at  the  composite  specimen  (or  structural)  level. 
The  two  analyses  are  conducted  in  tandem  as  damage  initiates  and  evolves  in 
composite  specimens  (structures).  Damage  is  described  at  the  micro-level  in  terms  of 
the  constituent  materials.  Damage  modes  for  a  particulate  composite  are  described  in 
terms  of  particle  cracking,  matrix  micro  cracking,  and  interface  debonding.  The 
interaction  between  the  two  levels  of  analyses  is  shown  in  (mg)  1  and  the  description 
of  each  analysis  is  given  below.  '  »  , 

■ 

Macro-Analysis  .  ,  l  .  „  a 


The  macro-analysis  is  conducted  at  the  composite  specimen  (structure)  level.  The 
specimen  considered  in  this  study  is  shown  in  Fig.  2.  It  is  a  square  plate  ^dth  a 
circular  hole  at  the  cjnter.  The  specimen  geometry  is  76.2  mm  x  76.2  mm  wftfr-6.35 
mm  thick.  The  raduses  of  the  hole  are  3.17pm  and  6.35  mm.  The  finite  element 
analysis  technique  is  used  for  the  macro-analysis.  Three-dimensional  solid  elements 
were  used  for  the  specimen.  Because  of  symmetry,  one^eighth  of  the  specimen  was 
modeled.  A  refined  mesh  was  used  around  the  notch  tip  to  capture  local  damage  and 
the  initial  crack. 
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Isoparametric  types  of  elements  are  used  with  the  Gauss  quadrature  rule.  These 
elements  require  material  properties  at  the  Gauss  quadrature  points.  As  a  result,  the 
interaction  between  the  micro-analysis  and  the  macro-level  finite  element  analysis  is 
performed  at  the  Gauss  quadrature  points.  The  micro-analysis  provides  the  Gauss 
quadrature  points  of  the  macro-analysis  with  effective  composite  materials  properties 
including  the  damage  state.  The  macro-analysis  is  undertaken  to  compute 
deformation  and  stress/strain  of  the  composite  specimen  (structure).  These  data  are 
transferred  to  the  micro-level  analysis. 


Micro- Analysis 


Because  the  micro-analysis  is  conducted  repeatedly  at  each  Gauss  quadrature  point  of 
finite  elements  as  damage  progresses,  the  analysis  should  be  computationally 
efficient.  To  this  end,  a  simplified  analytical  model,  or  a  unit-cell  model,  is  used  for 
the  micro-level  analysis. 

The  micro-analysis  has  two  major  functions.  One  is  to  compute  effective  material 
properties  of  the  composite  based  on  the  progressive  damage  at  the  constituent 
material  level.  These  data  are  used  for  the  subsequent  macro-analysis.  The  other 


function  is  to  decompose  the  macro-level  stress/strain  obtained  from  the  previous 
macro-analysis  into  micro-level  stress/strain.  The  micro-stress/strain  is  used  to 
determine  the  damage  state  at  the  constituent  material  level.  The  continuum  damage 
theory  is  utilized  to  describe  the  damage  state  as  well  as  the  degraded  material 
properties  of  each  constituent  material.  The  micro-analysis  utilizes  the  degraded 
material  properties  to  calculate  the  updated  material  properties  of  the  damaged 
composite. 

Continuum  Damage  Theory 

Damage  at  the  micro-level  is  described  using  a  continuum  damage  theory.  In  general, 
an  anisotropic  damage  theory  can  be  applied  to  the  constituent  materials  because  the 
damage  state  may  be  anisotropic  depending  on  the  given  conditions.  An  anisotropic 
damage  theory,  however,  requires  many  material  parameters^most  of  which  are  not 
easily  obtainable.  On  the  other  hand,  a  simpler  isotropic  damage  theory  may  be  used 
at  the  micro-level^ssuming  that  each  constituent  material  behavior  is  isotropic  at  the 
level  and  its  damage  state  is  also  isotropic.  These  assumptions  are  realistic.  Therefore, 
an  isotropic  damage  theory  is  used  in  this  study.  The  particulate  composite  under 
consideration  has  matrix  damage^  including  interface  debonding  but  not  particle 
cracking.  As  a  result,  the  discussion  below  applies  to  the  matrix  material. 

For  an  isotropic  damage  theory,  a  scalar  damage  parameter  is  used  to  indicate  the 
degree  of  damage  in  the  matrix  material.  The  damage  parameter  d  varies  from  0 
(denoting  no  damage)  to  dc  <  1  (denoting  damage  saturation).  When  damage 
saturates  in  the  material,  the  stress  in  the  material  reduces  to  zero  because  the  material 
cannot  support  the  external  load.  The  damaged  material  has  reduced  stiffness  as 
expressed  below: 

Eyki  =  EykiO-  ~  d)  (1) 

in  which  superscripts  V  and  ‘o’  denote  reduced  and  original  moduli,  respectively. 

The  damage  evolution  function  for  the  damage  parameter  is  described  based  on  an 
experimental  observation.  The  experimental  study  showed  that  damage  increases  as 
the  hydrostatic  stress  increases.  Therefore,  the  dilatational  strain  energy  was  an 
important  part  in  the  damage  process.  Previous  studies  used  the  total  strain  energy 
for  the  damage  evolution  function,  and  their  results  were  in  good  agreement  with  the 
experimental  data.  In  the  present  study,  the  damage  evolution  function  is  assumed  to 
be  a  function  of  dilatational  strain  energy  of  the  matrix  material.  Under  a  uniaxial 
load/vwhich  results  in  a  nearly  uniaxial  state  of  stress  near  the  notch  tip,  either  using 
the  total  strain  energy  or  the  dilatational  strain  energy  for  the  damage  function  may 
not  make  much  difference  because  both  functions  are  proportional  to  the  square  of  the 
uniaxial  stress  component.  However,  for  a  complex  loading,  the  total  strain  energy 
would  be  a  better  choice.  In  this  study,  the  rate  of  damage  increase,  d ,  is  assumed  to 


be  proportional  to  the  rate  of  dilatational  strain  energy  increase,  Ud ,  if  the  present 
dilatational  strain  energy  exceeds  the  threshold  value  U°d . 

d*Ud  if Ud>U°  (2) 

The  proportional  constant  and  the  threshold  dilatational  strain  energy  are  the  material 
properties  of  the  matrix.  The  damage  model  presented  here  includes  both  matrix 
cracking  and  interface  cracking  at  the  micro-level.  However,  if  the  interface  damage 
is  modeled  separately,  the  following  approach  can  be  used^but  it  was  not  used  for  the 
results  presented  in  this  paper.  One  of  the  reasons  the  interface  damage  was  included 
in  the  matrix  damage  is  that  interface  strength  data  is  not  available  with  any  associated 
failure  criterion.  However,  separate  modeling  of  matrix  damage  and  interface 
debonding  would  result  in  a  more  detailed  failure  process. 

For  particle/matrix  interface  damage,  it  is  assumed  that  the  damage  growth  at  the 
interface  is  related  to  the  change  of  the  volume  jSee' : 

d<Hl-<l>)dsu  (3) 

)5/here  <|>  is  the  interface  damage  parameter  and  suis  the  strain  dilatation.  This  equation 
applies  when  the  induced  strain  measure  exceeds  the  threshold  value  that  is  material 
dependent. 

Criterion  for  Initial  Crack  Length 

As  the  external  load  (uniform  displacement  in  this  case)  is  applied  to  the  specimen^ 
shown  in  Fig.  2,  there  is  a  stress/strain  rise  near  the  notch  tip  (see  curve  #1  in  Fig^Sa) 
before  damage  occurs.  The  curve  #1  in  Fig.  3  a  plots  the  characteristic,  induced  stfSss 
distribution  along  the  minimum  section  across  the  hole  of  the  specimen  before 
damage  initiation.  Damage  initiates  3 at  the  notch  tip  and  in  its  immediate 
neighborhood,  and(rt)grows  there  faster  than  at  other  locations  of  the  specimen.  Once 
damage  begins  to  accumulate  at  the  notch  tip  zone,  the  material  in  the  zone  becomes 
softens  described.  The  load-bearing  capability  of  the  softer  material  zone  reduces  by 
transferring  its  burden  to  the  neighboring  material  zone  with  less  damage  (i.e.  stiffer 
material).  Therefore,  the  induced  stress  distribution  along  the  minimum  section  is 
represented  in  curve  #2  of  Fig.^a  as  damage  grows  at  the  notch  tip  zone.  Further 
damage  growth  eventually  results^  in  a  lower  stress  at  the  zone  very  close  to  the  notch 
tip  than  at  its  neighboring  zone^sas  seen  in  curve  #3  of  Fig.  3  a.  Finally,  as  damage 
saturates  at  the  notch  tip,  the  induced  stress  vanishes  thereAas  shown  in  Fig.  3b.  No” 
load-carrying  capability  at  the  notch  tip  indicates  a  crack  formation  at  the  notch  tip  at 
the  given  load.  Then,  the  main  question  is  how  far  the  crack  will  propagate  to  be  the 
initial  crack  at  the  given  load. 
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In  order  to  determine  the  initial  crack  size,  the  unstable  (softening)  material  zone  is 
determined.  The  unstable  material  zone  means  the  induced  stress  decreases  even  if 
the  applied  strain  (displacement)  increases.  Once  the  crack  forms  at  the  notch  tip,  it  is 
considered  to  propagate  through  the  unstable  material  zone  until  it  meets  the  stable 
material  zone.  For  example,  Fig.3b  illustrates  both  unstable  and  stable  material  zones. 
In  summary,  the  length  ofAinitiail  crack  is  assumed  to  be  the  same  as  the  unstable 
(softening)  material  zone  when  damage  saturates  at  the  notch  tip  (i.e.  the  stress  at  the 
notch  tip  vanishes). 


Results  , 

Figure  4  shows  the  stress-stram/ curves  for  the  uncracked  specimen  tested  at  a 
constant  strain  rate  of  6.67  miny)  Also  included  is  the  predicted  stress-strain  curve 
based  on  the  finite  element  analysis.  It  is  seen  that  the  stress-strain  curves  obtained 
from  test  and  numerical  analysis  agre^sj  very  ^  well.  Similarly,  for  the  cracked 
specimen,  there  is  a  good  correlation  existj^etween  the  measured  and  the  predicted 
stress-strain  curves.  Figures  5  and  6  show  the  distributions  of  the  normal  stress  as  a 
function  of  the  normalized  distance  from  the  edge  of  the  hole.  According  to  these  two 
figures/\together  with  the  criterion  for  the  determination  of  the  initial  crack  length,  the 
predicted  initial  crack  lengths  are  0.0489  in.  and  0.0242  in.  for  the  small  (r  =  0.  125  ^  ^ 

in.)  and  the  large  size^fr  =  0.25  in.)  holef  respectively^^-irr^  The  measured  laveragddV 
values  of  the  initial  crack  length  are  0.0307  i%and  0.0526  in.  for  the  small  and  the  <*. 

large  sizej^-  holes,  respectively.  The  good  correlation  between  the  predicted  and  the  .  r-P 

measured  valuesAof  the  initial  crack  length^indicates  that  the  proposed  method  can  be  ^  ^ 
used  to  predict  the  initial  crack  length  with  good  accuracy. 


^  ...  a 

Conclusions  oe. 

*  .  _  • 

In  this  study,  a  micro-macro  approach  was  used  for  the  prediction  of  initiaL crack 
length  initiating  from  the  edge  of  the  hole  under  a  constant  strain  rate  condition.  The 
unstable  material  zone  caused  by  large  damage  was  defined  and  used  criterion 

for  the  initial  crack  size.  The  predicted  results  agreed  well  with  the  experimental 
measurements. 
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Figure  1.  Interaction  between  Micro-Analysis  and  Macro-analysis 
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Figure  2$  Finite  Element  Models 


Figure  3.  Stress  Distribution  gong  the  Minimum  Section  from  the  Notch  Tip  as  a  Junction  of 
damage:  (a)  Damage  Increases  from  Curve  1  to  Curve  3,  (b)  Stable  and  Unstable  Zones  when 

Damage  Saturates  at  the  Notch  Tip 
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Figure  4.  Stress-Strain  Curve 
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Figure  5.  Normal  Stress  Distribution  as  a  Function  of  the  Normalized 
Distance  from  the  Edge  of  the  Hole  (0.5  injfole  Diameter)^ 


